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Abstract 

The average charged track multipUcity and the normalised distribution of the scaled 
momentum, Xp, of charged final state hadrons are measured in deep-inelastic ep scattering 
at high in the Breit frame of reference. The analysis covers the range of photon vir- 
tuahty 100 < < 20 000 GeV^. Compared with previous results presented by HERA 
experiments this analysis has a significantly higher statistical precision and extends the 
phase space to higher and to the full range of Xp. The results are compared with e+e^ 
annihilation data and with various calculations based on perturbative QCD using different 
models of the hadronisation process. 
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1 Introduction 



The study of parton fragmentation and hadronisation processes provides valuable insights into 
the non-perturbative regime of Quantum Chromodynamics (QCD). These processes may be 
studied at HERA using inclusive charged particle production. In deep-inelastic ep scattering 
(DIS) the measurement of particle momentum spectra can be performed in the current hemi- 
sphere of the Breit frame [1], where the photon virtuality, Q, can be related to the momentum 
of the scattered parton. The charged hadron multiplicity and the distribution of their momenta 
scaled by (5/2 are the observables which are used in this analysis to study the fragmentation 
process in DIS. They can be directly compared with similar observables measured in one hemi- 
sphere of the hadronic final states in e+e" annihilation events, where particle momenta are 
scaled to half of the centre-of-mass energy E* /2. 

Previous comparison of DIS with several e+e^ experiments [2-5] have shown, in general, 
good agreement between these processes at high Q. At lower Q this agreement is observed 
to break down due to higher order QCD processes such as Boson Gluon Fusion (BGF) and 
Initial state Compton QCD (ICQCD). These processes occur as part of the hard interaction 
in ep scattering but not in e^e~ annihilation. They may lead to a relative depletion of the 
track multiplicity in the current region of DIS interactions as shown in [6]. Leading order 
matrix element Monte-Carlo programs, which use models of the parton cascade to describe 
QCD processes beyond leading order, have been shown to be able to describe the spectra in ep 
interactions down to low Q. 

Compared with the previous HI publication [3] this analysis utilises a ten times larger data 
sample with a better understanding of the experimental uncertainties. Results are obtained for 
a large range in Q (10 < Q < 100 GeV), which overlaps with several e~^e~ experiments 
including LEPl, and for the full range of the scaled charged hadron momenta. 



2 The HI detector 

A full description of the HI detector can be found elsewhere [7] and only the components most 
relevant for this analysis are briefly mentioned here. The origin of the HI coordinate system is 
the nominal ep interaction point, the direction of the proton beam defining the positive 2;-axis 
(forward region). 

The Liquid Argon (LAr) calorimeter measures the positions and energies of particles, in- 
cluding the scattered positron, over the polar angle range 4° < 9 < 154°. The calorimeter 
consists of an electromagnetic section with lead absorbers and a hadronic section with steel 
absorbers. The energy resolution for electrons in the electromagnetic section is a{E)/E = 
11.5%/VE [GeV] © 1% [8]. 

High events are triggered mainly using information from the LAr calorimeter. The trig- 
ger selects localised energy deposits in the electromagnetic section. For scattered positrons with 
energy above 11 GeV the trigger inefficiency is negligible as determined using independently 
triggered samples of events. There is no explicit track requirement in the trigger. 
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Charged particles are measured in the Central Tracking Detector (CTD) in the range 
20° < 9 < 165°. The CTD comprises two large cylindrical Central Jet Chambers (CJCs) 
arranged concentrically around the beam-line, complemented by a silicon vertex detector [9] 
covering the range 30° < 9 < 150°, two 2;-drift chambers and two multiwire proportional 
chambers for triggering purposes, all within a solenoidal magnetic field of strength 1.16 T. The 
transverse momentum resolution is a{p'r) /pr — 0.006 pt [GeV] © 0.02 [10] . In each event 
the tracks are used in a common fit procedure to determine the ep interaction vertex. 



3 Data Selection 

The data used in this analysis correspond to an integrated luminosity of 44 pb"^ and were taken 
by HI in the year 2000 when protons with an energy of 920 GeV collided with positrons with 
an energy of 27.5 GeV. 

Events are selected if the scattered positron is detected in the LAr calorimeter in the polar 
angular range 10° < < 150° and with energy greater than 11 GeV. The kinematic phase 
space, calculated using the scattered positron only, is defined by requiring the photon negative 
momentum Q'^ to be in the range 100 < < 20 000 GeV^ and the inelasticity y, defined as the 
fractional energy loss of the electron in the proton rest frame, to be in the range 0.05 < y < 0.6. 
The polar scattering angle for a massless parton, calculated from the positron kinematics in the 
quark-parton model (QPM) approximation \ is required to be in the range 30° < 9q iab < 150°. 
This ensures that the current region of the Breit frame remains in the central region of the 
detector where there is high acceptance and track reconstruction efficiency. It should be noted 
that the defined kinematic phase space can be applied simply to the theoretical models. 

Additional selections are made to reduce QED radiation effects and to suppress background 
events. The integrated hadronic final state is reconstructed from combined objects, built from 
calorimeter clusters and tracks, using an energy flow algorithm which ensures that no double 
counting of energy occurs. In order to minimise the correction due to QED radiation the value 
of y calculated from the hadronic final state using the Jacquet-Blondel method, yjh, and that 
calculated from the scattered positron, y^, are required to satisfy yjb — ye > —0.15 and {yjb — 
Vei/yjb > —0.75. The z coordinate of the event vertex is required to be within 35 cm of 
the nominal interaction point. This together with the rejection of events that have an event 
timing which does not match the HERA bunch crossing removes background from beam gas 
interactions. The difference between the total energy E and the longitudinal component of 
the total momentum Pz , calculated from the electron and the hadronic final state, E — Pz 
is required to be in the range 35 < — < 70 GeV in order to reduce background from 
photoproduction events. The event selection outlined above results in a data sample of about 
60,000 events. 

The reconstructed charged tracks in the selected events are used to study the fragmentation 
process. Only tracks that lie within the acceptance of the CTD which are fitted to the primary 

'The definition of the polar scattering angle is Oqjab = '^0'S~^( |^s(as-Q^)+4g^Q^ )' where E is the incoming 
positron beam energy, s is the ep centre of mass energy squared and x is the fraction of the proton momentum 
carried by the struck quark in the QPM. 
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vertex and have transverse momenta above 120 MeV are used in this analysis. In addition a 
variety of other track quality cuts are applied to remove badly measured tracks in a manner 
which can be accurately simulated. By using only tracks fitted to the event vertex the contribu- 
tion from the in-flight decays of A'°'s, A's, from photon conversions and from other secondary 
decays is minimised. 



4 Observables 

The Breit frame provides a kinematic region where the properties of the scattered quark can 
be studied with a well defined and relatively clean separation from the proton remnants. In the 
Breit frame of reference the virtual space-like photon has momentum Q but no energy. The 
photon direction defines the negative 2;'-axis and the current hemisphere. 

Within the QPM the photon collides head on with a (massless) quark of longitudinal mo- 
mentum Q/2. The struck quark thus scatters with an equal but opposite momentum into the 
current hemisphere while the proton remnants go into the opposite hemisphere. For the purpose 
of comparison the current region is taken to be the equivalent of one hemisphere of an e"^e~ 
annihilation. The energy scale, set by the virtual photon at (5/2, is taken to be equivalent to half 
of the e'^e" centre-of-mass energy E* /2. 

The boost to the Breit frame is defined using kinematics calculated from the properties of the 
scattered positron. Hadrons emerging from the interaction with negative longitudinal momenta 
in this frame are assigned to the current region and associated with the struck quark. 

Within this analysis the average charged multiplicity, <n>, is defined to be the average 
number of charged particles in the current region of the Breit frame per event. It is compared 
directly with half the average charged particle event multiplicity seen in e~^e~ annihilation. 

The scaled momentum variable Xp is defined to be Ph/iQ/ 2) where ph is the momentum of a 
charged track in the current region of the Breit frame. In e+e~ annihilation events the equivalent 
variable is ph/{E*/2). The inclusive, event normalised, charged particle scaled momentum 
distribution, D{xp, Q), is calculated as where in each Q range, is the total number of 

selected events and dn is the total number of charged tracks with scaled momentum Xp in the 
interval dxp. 



5 Phenomenology 

Fragmentation can be studied separately from the hard subprocess. The comparison of the 
fragmentation of the struck quark from the proton in DIS with that of a quark produced from 
e+e~ annihilation allows a test of quark fragmentation universality. Results from a number of 
different e~^e~ experiments [11, 12] at different centre of mass energies are available allowing 
comparison over the full Q range of this analysis. The contribution from weakly decaying short 
lived particles (e.g. and A) is subtracted from the e~^e~ results to be consistent with our 
charged particle selection. This contribution is about 8% and is estimated from the ep Monte 
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Carlo program DJANGO [13] which provides in general a good description of strange particle 
production [14]. 

The data presented here are used to test predictions of different models of the parton cascade 
and hadronisation processes, implemented in various leading order matrix element Monte Carlo 
programs, which have been tuned to describe e~^e~ results. 

The Parton Shower model (PS) [15], implemented in the RAPGAP [16] Monte Carlo pro- 
gram, describes the fragmentation process as the splitting of a parent parton into two daughters 
(q ^ qg, g ^ gg, g qq)- The splitting continues, giving rise to a parton shower. The 
evolution of the parton shower is based on leading log Q"^ DGLAP [17] splitting functions. The 
transverse momentum, kr, of subsequently emitted partons is highly ordered. Gluon coher- 
ence, which suppresses the emission of soft gluons at wide angles, is approximately modelled 
by imposing angular-ordering. 

In the Soft Colour Interaction model (PS-i-SCI) [18] soft gluons are exchanged between the 
partons produced in the parton shower and the proton remnant. Soft Colour interactions are 
simulated using the implementation available in the LEPTO [19] Monte Carlo program. 

In the Colour Dipole Model (CDM) [20], dipoles are created between coloured partons. 
Gluon emission is treated as radiation from these dipoles. New dipoles are formed with the 
emitted gluons from which further radiation is possible. The radiation pattern of the dipoles in- 
cludes interference effects, thus modelling gluon coherence. The of emitted partons are only 
weakly ordered, producing a picture similar to the BFKL treatment of parton evolution [21]. 
ARIADNE [22] provides an implementation of the colour dipole model and is used in the 
DJANGO [13] Monte Carlo program. 

The RAPGAP, LEPTO and DJANGO Monte Carlo programs use the Lund string model of 
hadronisation [23] which is based on the dynamics of a relativistic string, or gluonic "flux tube", 
stretched between coloured partons. As the partons move apart they lose kinetic energy to the 
string creating pairs which form new string pieces. This process is iterated until the available 
energy is used up. The residual string fragments are combined into mesons and baryons. 

The HERWIG Monte Carlo [24] program uses the parton shower model to describe the 
fragmentation process but incorporates the cluster model of hadronisation [25]. In the cluster 
model partons are generated in a perturbative shower. The cascade is stopped at a given cut-off, 
related to the minimum transverse momentum of the emitted partons. The remaining gluons are 
then split into light qq pairs. Coloured objects that are close to one another are combined into 
colourless clusters which decay isotropically in the rest frame of the original cluster into known 
resonances. 

Next-to-leading order (NLO) perturbative QCD (pQCD) predictions based on the CYCLOPS 
program [26]are also available. In CYCLOPS the DIS cross-section is factorised into three 
parts: a proton parton-density function (PDF), the full NLO matrix element (ME), and a par- 
tonic fragmentation function (EE). The PDE parameterisation CTEQ6M [28] is taken as default 
and the results are cross-checked using MRST2001 [29]. The factorisation and renormalisation 
scales are chosen as Q. The uncertainty arising from the scale choice is estimated by increas- 
ing and decreasing the scale by a factor of two. Three different parameterisations of the EE 
are used which are obtained from NLO fits to e+e" data: KKP [30], KRETZER [31], and 
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AKK [32]. Each parameterisation uses the appropriate quark mix (u,d,s,c,b) for ep interactions. 
This mix is not flavour democratic unlike the case of e+e~ annihilations. The predictions are 
for light charged hadron production (tt^, K^, and 'p') as measured in this analysis. The value 
of Ai|I = 266 MeV is used for the PDF, ME and EE. 

In order to avoid infrared singularities in the NLO calculation, the hadrons should be pro- 
duced close in rapidity to the parent parton. Therefore the rapidity of the hadron in the Breit 
frame is required to be greater than unity [27]. Within this analysis this corresponds to a safe 
limit of Xp > 0.1. As a consequence the data are only compared to the NLO QCD predictions 
for Xp > 0.1. 

6 Data Correction 

The data are corrected for detector acceptance, efficiency and resolution effects using Monte 
Carlo event samples generated with the RAPGAP and DJANGO programs. All generated 
events are passed through the full GEANT [33] based simulation of the HI apparatus and are 
reconstructed and analysed using the same programs as for the data. These Monte Carlo event 
samples give a good description of the data. The residual contribution of charged particles from 
the weak decay of neutral particles (e.g. and A's) is subtracted from the data as part of the 
correction procedure. The effects of QED radiation are corrected for using the HERACLES [13] 
program incorporated within the above Monte Carlos. The total correction factor is calculated 
from the ratio of the number of entries in each bin at hadron level to that at detector level. The 
bin sizes are chosen to give high acceptance and purity^, typically above 60%, with a minimum 
of 40%. The total correction factor applied to the uncorrected data points is ~1.1 for <n> and 
typically less than 1.2 for D(xp, Q). In general the uncertainty in the boost to the Breit frame 
dominates the resolution in Xp. The correction associated with the tracking dominates in the 
highest Q'^ region where there is a somewhat reduced acceptance for the current region of the 
Breit frame within the CTD. 



7 Systematic Uncertainties 

The following sources of systematic errors are considered for all measured quantities. Also 
presented are the resulting typical fractional error on <n> and, where different, D(xp, Q). 

• The positron energy scale uncertainty is 0.7 — 3 % depending on the position of the 
detected positron in the LAr calorimeter. This uncertainty affects both the phase space 
and boost calculation. Its effect is studied by repeating the analysis while varying the 
positron energy scale in the simulation used to correct the data. This gives an error on 
<n> of about 2 % independent of Q. The resulting uncertainty on D{xp, Q) is again 
independent of Q but varies with Xp from 0.5% {Xp ~ 0.1) to 7% {Xp ~ 1.0). The 
scattered positron angular resolution leads to a systematic error of about 1%. 

^The acceptance (purity) is defined as the ratio of the number of charged hadrons generated and reconstructed 
in the bin to the total number of charged hadrons generated (reconstructed) in that bin. 
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• The systematic error associated with the track reconstruction efficiency is estimated to 
be 2.5%. This is applied as an independent uncorrelated error on every point and is 
assumed to be constant over the pt and 6 range of the tracks. In the highest Q interval 
an additional uncertainty due to problematic track reconstruction in dense jets with high 
track multiplicities leads to an extra error of 2% on the average charged multiplicity and 
5%onD(xp,g). 

• The hadronic energy scale uncertainty is taken to be 4%. Only the calculation of E — Pz 
and iijb, used in the event selection, rely on the hadronic energy measurement. Varying 
the hadronic energy scale in the simulation of the samples used to correct the data gives 
an error of about 1%. 

• The uncertainty in the correction factor arising from using different Monte Carlo models 
in the correction procedure, taken as the full difference between correcting the data with 
RAPGAP or DJANGO, resuks in a typical error of 1% on <n> and 1.5% on D{xp, Q). 

Apart from the error from the scattered positron energy, which is quoted separately, the 
individual effects of the above experimental uncertainties are combined in quadrature. The total 
systematic uncertainty is dominated by the systematics attributed to the tracking efficiency. In 
the defined kinematic region errors arising from non ep background are negligible. 



8 Results 

The measurements of the averaged charged particle multiplicity, <n>, and the scaled momen- 
tum distribution, D{xp,Q), are presented. The residual and A contribution is subtracted 
from all results. The data are listed in tables 2 and 3 and shown in figures 1 to 5 at the average 
Q values given in table 1. 

8.1 Average Charged Multiplicity 

In figure la) the measurements of <n> are compared with different parameterisations of the 
average charged track multiplicity per event seen in e~^e~ annihilation taken from [11] and with 
results from the ZEUS experiment [5]. 

The ZEUS results are in agreement, within the errors, with the HI data. For most of the 
Q range the HI data are in good agreement with the parameterisation of the e+e" data. In 
the lowest Q interval the data are slightly below the parameterisation of the e"'"e~ data. In the 
highest two Q intervals the measurements are clearly below the e~^e~ parameterisation. 

In figure lb) a comparison is made with different models of the hadronisation and parton 
cascade processes implemented in leading order matrix element Monte Carlo programs. The 
models show good agreement with our data except for the one with soft colour interactions 
which significantly overestimates the multiplicity. 
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Extensive tests have verified that the difference between the HI data and the e~^e" param- 
eterisation is not due to any sensitivity to phase space selection. At a given Q, <n> varies by 
less than 2% over the x, y or Og^iab ranges. Moreover, in the vicinity of the Z° resonance the 
admixture of heavy and light flavoured quarks in ep and e+e^, especially with respect to the 
amount of b quarks, gives an overall excess in the e~^e~ multiplicity of ~ 3% which can not 
explain the observed difference. A possible explanation of this difference is the complexity of 
the DIS scattering, with extra colour connections between the scattered quark and the proton 
remnants, which complicates the simple analogy with e+e" data. The Monte Carlo models, 
which include some of the additional processes reflecting the complexity of the DIS interaction, 
are compatible with the data apart from the SCI model. 

8.2 Scaled Momenta Spectra 

In figure 2 the normalised distribution of the scaled momentum, j^^, is shown as a function 
of Q for nine different intervals of Xp. The results are compared to e'^e~ annihilation data [12]. 

At low Xp moving from low to high Q there is a dramatic increase in the number of hadrons 
while at high Xp the number of hadrons decreases with Q, i.e. the spectra becomes softer as Q 
increases. This well known observation contradicts the scaling hypothesis of the fragmentation 
functions, which would imply an independence of Q of the scaled momenta of hadrons. The ep 
data show similar behaviour to the e+e~ data providing a rough demonstration of fragmentation 
universality. In certain phase space regions some discrepancies are visible. It has been observed 
previously [3-5] that for Q < 10 GeV there is a discrepancy between ep and e^e~ data. This 
is understood in terms of higher order QCD processes depleting the current region as described 
in the introduction. In this analysis a similar difference is seen at low Q (Q ~ 15 GeV) and 
intermediate Xp (0.05 < Xp < 0.4). This is reflected by the observation of a slightly reduced 
<n> for the lowest Q interval, and suggests that these higher order QCD processes have an 
influence in this range. At high Q {Q > 60 GeV) and small values of Xp (0.02 < Xp < 0.2) 
there are significantly less tracks observed in the ep data compared to e'^e~ annihilation. This 
corresponds to the observed <n> values below the expectation in this Q range as described in 
the previous section. 

In order to investigate the fragmentation process in greater detail it is necessary to go beyond 
a simple comparison with e+e~ data and compare with models that account for the processes 
specific to ep scattering. In figure 3 the data are compared with Monte Carlos model predictions 
that implement different models to describe the parton cascade and the hadronisation process. 

It can be seen from figure 3 that both CDM and the PS model provide an acceptable descrip- 
tion of the data, with the CDM model predicting a slightly harder spectrum. Both models tend 
to overestimate the multiplicity at higher Q. The SCI model predicts a much softer spectrum 
than the other two models and is disfavoured by the data. This could be due to the additional 
gluon interactions in the SCI model which soften the spectra of partons produced by the parton 
shower. 

The data clearly show a preference for predictions where the Lund string model of hadroni- 
sation is used (PS and CDM). HERWIG (cluster hadronisation) predicts a spectrum that is too 



10 



hard compared to the data. At high Xp HERWIG even fails to reproduce the observed scaling 
violation and predicts a flat spectrum in Q. 

In figure 4 the data are compared with the predictions obtained from the NLO Monte Carlo 
program CYCLOPS [26] for three different parameterisations of the FF in the infrared safe 
region as defined in section 5. The uncertainties associated with the change in scale or PDF 
are significantly smaller than the differences between these parameterisations. The different 
fragmentation function parameterisations give different results but it is evident that none of them 
can describe the scaling violations seen in the data. This prevents a reliable extraction of either 
tts or of the fragmentation function from the data measurements using this NLO calculation. It 
is interesting to note that the KKP and KRETZER fits use different assumptions for the light 
quark flavour contributions to the fragmentation function while the AKK fit uses recent e^e^ 
data which include light quark tagging probabilities to constrain the strangeness contribution. 
The AKK and KRETZER parameterisations are seen to agree quite well with each other. 

A summary of the results for the scaled momentum spectra is presented in figure 5a) as a 
function of Xp for different Q intervals, where each Q interval has been scaled by an additional 
factor ten for improved visual display, and are compared with the PS Monte Carlo prediction. In 
figure 5b) the results are presented as a function of Q for different Xp intervals and are compared 
with the e+e^ annihilation data and the PS Monte Carlo prediction. 



9 Conclusions 

The average charged multiplicity, <n>, and the event normalised scaled momentum distribu- 
tion, D{xp, Q), of charged hadrons have been measured in ep collisions at high in the Breit 
frame of reference and compared with e~^e~ data and a variety of models. 

The results broadly support the concept of quark fragmentation universality in ep collisions 
and e"'"e" annihilation. A small multiplicity depletion compared to e~^e~ is observed at low Q 
which can be attributed to higher order QCD processes occurring as part of the hard interaction 
in ep scattering but not in e"^e" annihilation. At high Q a large depletion is observed. 

The best description of the data by leading order matrix element Monte Carlo programs 
is given by models that use the string model of hadronisation and do not include soft colour 
interactions. In the low and high Q regions, where the comparison to e"^e^ is poor, the Monte 
Carlo models are able to provide a better description of the data. 

The results are compared with NLO QCD calculations as implemented in the CYCLOPS 
program. All three parameterisations of the fragmentation functions used in this program fail to 
describe the scaling violations seen in the data. 
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(GeV^) 


<g> GeV 


oQ GeV 


<x> 


ox 


100 < g2 < 175 


12.3 


0.1 


0.00370 


0.00004 


175 < g2 < 250 


14.5 


0.1 


0.00952 


0.00007 


250 <Q'^ < 450 


18.0 


0.1 


0.1559 


0.0001 


450 < g2 < 1000 


25.0 


0.3 


0.0254 


0.0003 


1000 < g2 < 2000 


36.6 


0.8 


0.044 


0.001 


2000 < g2 < 8000 


58.5 


2.1 


0.087 


0.003 


8000 < g2 < 20000 


102.0 


17.0 


0.20 


0.03 



Table 1 : Average Q and x values and their statistical errors for the selected events in the 
intervals used in this analysis. 



g2 (GeV2) 


<n> 


^stat 


5tot 


^ scale 


100 < g2 < 175 


3.39 


0.02 


0.09 


0.05 


175 < g2 < 250 


3.89 


0.01 


0.10 


0.08 


250 < g2 < 450 


4.41 


0.01 


0.12 


0.07 


450 < g2 < 1000 


5.21 


0.03 


0.13 


0.07 


1000 < g2 < 2000 


6.22 


0.06 


0.17 


0.07 


2000 < g2 < 8000 


7.37 


0.12 


0.22 


0.09 


8000 <Q'' < 20000 


8.11 


0.47 


0.55 


0.24 



Table 2: Average charged particle multiplicity <n> as a function of shown with the statis- 
tical error {5stat), the total error including statistical and systematic errors added in quadrature 
{5tot ), and the correlated error coming from the electron energy scale uncertainty {5scaie) which 
is not included in the total error. 
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Table 3: The measured normalised distribution of the scaled momentum j^dn/dxp as a function 
of for different Xp intervals shown with the statistical error (^stat)^ the total error including 
statistical and systematic errors added in quadrature {5tot and the correlated error coming from 
the electron energy scale uncertainty {5 scale) which is shown as two numbers (+/— ) and is not 
included in the total error. 
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Figure 1: The average charged multiplicity as a function of Q. For each measurement the 
statistical error is shown by the inner error bar while the outer error bar represents the statistical 
and systematic errors added in quadrature. In addition there is a further correlated error of 2% 
coming from the electron energy scale uncertainty (not shown). The data are displayed at the 
average value of Q, the horizontal error bars represent the statistical errors in table 1. The data 
are compared with: a) parameterisations of data from e'^e' experiments [12] (taking Q = E*) 
and with results reported by the ZEUS experiment, and b) predictions from different models of 
the hadronisation and parton cascade processes implemented in leading order matrix element 
Monte Carlo programs as described in the text. 



16 



< Xp < 0.02 



|"200 
2 150 
100 
50 



• H1 Data ^ 

★ DELPHI '■^ 
T TASSO 

▲ MARKII « 

■ AMY 



10 



10^ 

Q,E* (GeV) 



0.02 < Xp < 0.05 



X 
7S 




10 



10' 

Q,E* (GeV) 



35 



0.05 < Xp < 0.1 



X 

C 30 ■ 
7S 



10 



25 h f 



20 ■ 
15 ■ 



10 



10^ 

Q,E* (GeV) 



0.1 < Xp < 0.2 



X 




10 



10" 

Q,E* (GeV) 



X 

■o 



0.2 < Xp < 0.3 




10 



10' 

Q,E* (GeV) 



0.3 < Xp < 0.4 



X 
73 




1.5 ■ 



10 



10' 

Q,E* (GeV) 



0.4 < Xp < 0.5 



X 

■o 



0.8 



0.6 



0.4 



i T 




i 


> 


■ i 



10 



10^ 

Q,E* (GeV) 



a. 
X 
73 

0.4 

73 



0.3 



0.5 < Xp < 0.7 



0.2 











T W -r J 

A It 


i 



10 



10^ 

Q,E* (GeV) 



0.7<Xp<1.0 




0.02 ■ 



10 



10^ 

Q,E* (GeV) 



Figure 2: The measured normalised distributions of the scaled momentum, j^dn/dxp, as a 
function of Q for nine different Xp regions. The statistical error is shown by the inner error bar 
and the outer error bar represents the statistical and systematic error added in quadrature. In 
addition there is a further correlated error of ~ 0.5 — 7% (increasing with Xp) coming from the 
electron energy scale uncertainty (not shown). The data are displayed at the average value of Q, 
the horizontal error bars represent the statistical errors given in table 1 . Data are compared to 
results from various e^e^ experiments (taking Q = E*). Note the suppressed zeros and large 
change in scale of the vertical axis moving to higher values of Xp. 
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Figure 3: The measured normalised distributions of the scaled momentum, j^dn/dxp, as a 
function of Q for nine different Xp regions. The statistical error is shown by the inner error bar 
and the outer error bar represents the statistical and systematic error added in quadrature. In 
addition there is a further correlated error of ~ 0.5 — 7% (increasing with Xp) coming from the 
electron energy scale uncertainty (not shown). The data are displayed at the average value of 
Q, the horizontal error bars represent the statistical errors given in table 1 . The data are com- 
pared to predictions from different models of the parton cascade and hadronisation processes 
implemented in leading order matrix element Monte Carlo programs as described in the text. 
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Figure 4: The measured normalised distributions of the scaled momentum, j^dn/dxp, as a 
function of Q for six different Xp regions where there exist infra red safe NLO QCD predic- 
tions. The statistical error is shown by the inner error bar and the outer error bar represents 
the statistical and systematic error added in quadrature. In addition there is a further correlated 
error of ~ 0.5 — 7% (increasing with Xp) coming from the electron energy scale uncertainty (not 
shown). The data are displayed at the average value of Q, the horizontal error bars represent 
the statistical errors given in table 1 . The data are compared to NLO QCD CYCLOPS predic- 
tions for Q < 60 GeV using three different fragmentation functions: KKP (dot-dashed line), 
AKK (solid), and KRETZER (dashed). The typical scale uncertainties for the AKK predictions 
are shown as a shaded band. As standard the CTEQ6.1 PDF is used, the effect of using the 
MRST2001 PDF for the AKK predictions are also shown (dotted). 
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Figure 5: The measured normalised distributions of the scaled momentum j^dn/dxp'. a) as a 
function of Xp for the different Q intervals compared with the PS Monte Carlo prediction. Each 
Q interval, apart from the lowest, has been scaled by an additional factor of ten; and b) as a 
function of Q for the different Xp intervals compared with the e^e~ annihilation data and the 
PS Monte Carlo prediction. 
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